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Abstract 

We review recent progress in the mathematical theory of quantum 
disordered systems: the Anderson transition (joint work with Domin- 
gos Marchetti) , the (quantum and classical) Edwards-anderson spin-glass 
model and return to equilibrium for a class of spin glass models, which 
includes the EA model initially in a very large transverse magnetic field. 

1 - Introduction 

In recent years there has been a significant progress in the mathematical 
theory of (quantum) disordered systems. Our purpose here is to present the 
main ideas (without proofs), with a clear discussion of their conceptual and 
physical relevance, as well as a brief comparison with the recent, analogous 
results in the literature. Some new remarks and results, which clarify some 
important conceptual points, arc also included. 

We shall be interested in properties of disordered systems at low tempera- 
tures, i.e., near the ground state (absolute temperature T = 0); in the critical 
region of spin glasses, there exist the spectacular recent rigorous results on the 
mean field theory (see the review by F. Guerra |Gue06j ). and various results for 
other disordered systems, including the well understood high temperature phase 
in spin glasses, are discussed in the comprehensive book by A. Bovier |Bov06| . 

Our restriction to very low temperatures implies, of course, that we shall be 
dealing exclusively with quantum systems. In particular, the Ising model, when 
it appears, must be regarded as the anisotropic limit of quantum (spin) systems: 
the fact that this is not only mathematically so is emphatically demonstrated 
by the well-known fact that the critical exponents of the Ising model in three 
dimensions (see, e.g., jZJ79j ) are surprisingly close to those measured in real 
magnetic systems, precisely because most of the latter are highly anisotropic. 

Three important issues appearing in the above-mentioned context are local- 
ization, arising in connection with the Anderson transition, which we discuss in 



section 2, reporting on joint work with Domingos Marchetti ( |MW12a] .[MW12b ), 
frustration, which appears in short-range spin glasses such as the Edwards- 
Anderson (EA) model, discussed in section 3, and the (thermodynamic version 
of) instability of the first and second kind, which relates to the return to 
equilibrium for special initial states and probability distributions in a class of 
models (including the EA model) in section 4. 
2 - The Anderson transition 

The breakdown of translation invariance (or, more generally, Galilean invari- 
ance in many-body systems with short-range forces) leads to the existence of 
crystals, and the corresponding Goldstone excitations are phonons |Swi67j . On 
the other hand, the discrete translation invariance subgroup of the crystal is also 
frequently broken by e.g. impurities. This fact brings about a number of im- 
portant new conceptual issues, and is usually modelled by the introduction of a 
random local potential in a tight-binding model (for the latter see |Fe"y65 1 ) . The 
associated physical picture consists of lightly or heavily doped semiconductors 
(e.g. Si doped with a neighboring element which contributes excess electrons, 
for instance P). This is the Anderson model ([And58 , And72 ), described by 
the Hamiltonian 

H w = A + V u 

(1) 

on l 2 (Z d ) where A is the (centered discrete Laplacian) 

( Au )n = Un ' 
n'\\ri— n'\ — 1 

(2) 

plus a perturbation by a random potential 

(V u u) n = V^u n 

(3) 

where {V^} n€Z d is a family of independent, identically distributed random vari- 
ables (i.i.d.r.v.) on a probability space (f2,#, fi), with a common distribution 
F(x) — [i ({uj : V£ < x}); V£ is assumed to depend linearly on a quantity v > 
(see the forthcoming (5b)), which is the disorder parameter, also called cou- 
pling constant. The spectrum of H u is, by the ergodic theorem, almost surely 
a nonrandom set a(H u ) = [—2d, 2d] + supp dF . Anderson conjectured that 
there exists a critical coupling constant < v c < oo such that for v > v c the 
spectral measure of H is pure point (p.p) for /^-almost every w, while, for v < v c 
the spectral measure of H u contains two components, separated by so called 
mobility edge E ± : if E E [E~ , E + ] the spectrum of H 1 ^ is pure absolutely con- 
tinuous (a.c); in the complementary set a(H u ')\[E^ , E + ], has pure point 
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spectrum - leading to the new important phenomenon of localization, first 
proved for d — 1 in GMP77,, and for d > 3 for the first time in by |JF85j . based 
on previous seminal work by |FS83| . for fixed energy and large disorder, or for 
fixed disorder and large energy (the latter meaning near the edges of the band, 
i.e, the boundary of a(H u ))\ see also |AM 93 for a much simpler proof by an 
entirely different method, and references given there. 

As is well known, resonance is a phenomenon of cooperation between two or 
more elements, while randomness is based on independence - noncooperation 
- between elements. As remarked by Howland }How92] . it should not prove sur- 
prising, therefore, that resonance in a system may be removed by randomizing 
parameters in it, preventing, in the Anderson model, the cooperation necessary 
for travelling waves (tunneling). What is, however, surprising, is the (conjec- 
tured, see [And72 ) sharpness of the transition, which means that for v — v c 
any slight (positive) variation of the disorder parameter destroys the tunneling. 

This tunneling instability was studied by Jona-Lasinio, Martinelli and 
Scoppola I..MSS 1 1 and Simon |Sim85j . who showed that tunneling is very sensi- 
tive to minimal changes in a double-well potential, even those localized very far 
away from the minima ( a phenomenon called "flea on the elephant" by Simon). 
Tunneling instability may be of dynamical nature (see, e.g., |GS01| . |WC98| ) . 
being also at the root of the existence of the chiral superselection rules induced 
by the environment which account for the shape of molecules such as Ammonia 
NH3, see the review by Arthur S. Wightman, "Superselection sectors: old and 
new", |Wig95| . We shall see that for a lattice model such as (1), the idea that 
only a very "slight" perturbation of the Hamiltonian has a drastic effect on 
the tunneling acquires even a new dimension in the sparse situation described 
below. 

The proof of the existence of a mobility edge in the Anderson model re- 
mains as one of the major open problems in mathematical physics. Given the 
difficulties in proving this for the model (l)-(3), one might be led to study the 
limit v of H, for which the spectrum is pure a. a. We shall instead follow a 
different approach to the Anderson conjecture suggested by Molchanov ([Mol98 
[Mol99 , [MV98J: the limit of zero concentration, i.e., taking V u in H such that 

v?=xy(»-oi) > 

i 

with elementary potential (" bump") (p u : Z d — > K satisfying a uniform inte- 
grability condition 



1^)1 < 



for some e > and < Cq < 00 and 



C\ 







l + \z\ d+E 



hm #{*:N<^}= 
R-yoa R d 



(4) 



3 



Due to condition of zero concentration, potentials such as V are called sparse 
and have been intensively studied in recent years since the seminal work by 
Pearson [Pea78j in dimension d = 1 and by Molchanov Mol98], Mol99 in the 
multidimensional case. For d > 2 the interaction between bumps is weak while 
for d — 1 the phase of the wave after propagation between distant bumps be- 
comes " stochastic" . This is the right moment to introduce our one-dimensional 
model. 

We consider (infinite) Jacobi matrices 

{J u) n = U n+ i + U„_! , 

with the perturbation potential 

{V u u) n = v>„ (5a) 

with u = (u n ) n >o e h(Z+) ,and 

u _ j v if n = e A, 
71 1 if otherwise , 

(56) 

for v G (0, 1). A = {a" }j>i denotes a random set of natural numbers a" = 
dj + ujj , where aj satisfy the " sparseness condition" : 

a i - a i-i = P 3 ) i = 2,3, ... 

with a\ + 1 = p > 2 and u>j, j > 1, are independent r.v. on a probability 
space (O, B, fi) uniformly distributed on a set Aj = {— j, . . . , j}. These variables 
introduce an uncertainty in the positions of those points for which v n ^ 0: such 
models are called Poisson models (see |Jit07| and references given there). 
Note that the support of the ujj only grows linearly with the suffix j. We write 
J" = Jo + V u , and denote by the operator associated to J u on the Hilbert 
space H. of the square integrable sequences u = {u n ) n> _ 1 which satisfy a <p-h.c. 
at -1: 



U—i cos (j) — uq sin = 

The essential spectrum of equ als [—2,2]: it will be represented as A = 
2 cos a with a £ [0, 7r). Zlatos [Zla04] proved that this model exhibits a sharp 
transition from s.c. to p.p. spectrum. This was shown independently in 
QdCMWll| . [MW12hp : 

Theorem 1 Let be as above. Let 

/= |ag [-2,2]\2costtQ : i(0 - 1)(4 - A 2 ) > l| 

with v € (0, 1) and /3 € N, /? > 2 . Then, for almost all uj with respect to the 
uniform product measure on A = X^Li {— j, ■ ■ ■ ,j}, 
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(a) there exists a set A\ of Lebesgue measure zero such that the spectrum 
restricted to the set I\A\ is purely singular continuous, 

(b) the spectrum of Jp^ is dense pure point when restricted to 

I c = ([—2, 2]\2cos7rQ) \I for almost every <f> G (0, 7r), where <f> characterizes the 
boundary condition. Thus it is purely p.p. in this interval. 



In the above figure, we see the s.c. (light gray) and p.p. (dark grey) spectra 
separated by the "mobility edges" X ± = ±2-^(1 - with ^ = 1, 3038 
Remark 1 

The basic emphasis on one-dimensional (sparse) models has an important 
technical reason: the profound approach of Gilbert and Pearson connecting 
the space asymptotics of eigenfunctions of the restrictions of a large class of 
Sturm-Liouville operators to finite intervals to the spectral theory of these same 
operators in infinite space through the concept of subordinacy [GP87] is only 
available in one dimension. Using the important transfer matrix version of this 
theory due to Last and Simon [LS99 , the surprise is that, even in a regime of 
strong sparsity, a spectral transition from s.c. to p.p. spectrum (first shown 
by Zlatos [Zla04] ) may be proved: this is theorem 1. Our proof (theorem 4.7 
of [MW1 2b ) differs from that in |Zla04] by the use of the (optimal) metrical 
version of Weyl's theorem on uniform distribution |KN74] (in this case, of the 
so-called Priifer angles, see |KLS98j ) due to Davenport, Erdos and LeVeque 
[DEL63] . The latter permits an explicit characterization of the exceptional set 
in which the spectrum may not be of p.p. type (this exceptional set arises in 
connection with the concept of essential support of a measure [GP871). which 
reveals it to be also of p.p. type. As a result, the spectrum at the edges of the 
band is proved to be purely p.p.. 

The robustness of the transition depicted in theorem 1 follows from [dRJLS96] , 
because the Hausdorff dimension of the s.c. spectrum may be seen to be nonzero 
f [Zla04] . jdCMWlOj L 

Finally, the difference Laplacean (2) leads to tunneling and bands |Fey65| , 
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with a pure a.c. specrrum. The perturbation (5) which leads to the emergence 
of p.p. spectrum is very slight on two different counts: 1) the transition arises 
at a sharp (critical) value v c of the disorder parameter v, 2) at the same time, 
we are in a regime of strong sparsity! In this sense, the analogy to the "flea 
on the elephant" |Sim85) in a double-well potential is here much sharper than 
in the originally conjectured Anderson transition, and the transition depicted 
in theorem 1, as well as in the higher dimensional version in theorem 2, should 
indeed be viewed as a concrete manifestation of tunneling instability. 

Unfortunately the s.c. spectrum does not posess either the dynamic or the 
perturbation theoretic properties (see |S W86J , How80 for the latter) which are 
commonly associated with the physical picture of delocalized states. Regarding 
the dynamical properties, for instance, the sojourn time of a particle with 
initial state $ in a compact region S is defined as 

/oo 
\\P s exp(~itH)^\\ 2 dt 
-oo 

where Ps denotes the projector on S. If ^ e Use, it follows from a theorem by 
KB. Sinha |Sin77j that a S exists with |5| < oo and J(5;*) = oo. See also 
|MW12aj . [MW12bj for further discussion. 

We try therefore to attain higher dimensions. 

Multidimensional version 

Consider the Kronecker sum J% as an operator on H^H: 
4 2) := Jf ®I + 6I® jf 

(6) 

where w 1 = (wj) . >x and cj 2 = (u^) are two independent sequences of inde- 
pendent random variables defined in (£l,B, u), as before (we omit to 1 and lu 2 in 
the l.h.s above for brevity). Above, the parameter 9 £ [0, 1] is included to avoid 

(2) 

resonances . We ask for properties of Jg (e.g. the spectral type) which hold for 
typical configurations, i.e., a.e. (u; 1 ,^ 2 ,^) with respect to v x v x I where I is 

the Lebesgue measure in [0, 1]. Jg is a special two-dimensional analog of J^; if 
the latter was replaced by — A + V on J 2 (M, dx) where A = d 2 /dx 2 is the second 
derivative operator, and V a multiplicative operator ViJ)(x) = V(x)ip(x) (po- 
tential), the sum above would correspond to {—d 2 /dx\ + Vi) + (—d 2 /dx\ + V2) 
on L 2 (M. 2 , dx\dx2), i.e., the "separable case" in two dimensions. Accordingly, 

(n) 

we shall also refer to Jg , n — 2, 3, . . ., as the separable case in n dimensions. 
Our approach is to look at the quantity 
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($,e"^ 2) *) =f\t)f{6t) 
where 

f(s) = f sc (s) + r pp (s), i = l,2 

with/* G ( S ) = [ e- lsX d^(\) 



Above 



$ = (Vl+Pl) O (<^2+P2) , 
* = (V'l+Xl) ® (V>2+X 2 ) , 

with <fi,ipi G ^sc, Pi,Xi € "Hpp and y+p denotes the direct sum of two vectors 
ip,pEH. 

We shall use the following folklore proposition (see, e.g., theorem 3.2 of 
I.MW 12bl or [5m77p : 
Proposition 1 

Let fi be a measure on the space M(K) of all finite regular Borel measures 
on R. If the Fourier-Stieltjes transform of \i 

R3t^ fi{t) = J e- ux d^{\) 

belongs to L 2 (M.,dt), then [i is absolutely continuous with respect to Lebesgue 
measure. 

The time-like decay of the Fourier-Stieltjes (F.S.) transform of the spectral 
measure is dictated by the Hausdorff dimension of the spectral measure: 
Let S be a subset of R, a e [0, 1], and S > 0. Define 

oo 

^(S) = inr{£|Q|* 

IScU^ljCj with \C S \ < 6} 
where |C| denotes the Lebesgue measure (length) of C, and 
h a (S) = lim h%{S) = sup h?(S) 

8>0 

We call h a a - dimensional Hausdorff measure on R. h 1 agrees with Lebesgue 
measure, and h° is the counting measure, so that {/i"|0 < a < 1} is a family 
which interpolates continuously between the counting measure and Lebesgue 
measure. 
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Definition 1 A Borel meaure on R is uniformly alpha-Hoelder contin- 
uous (UaH) if there exists C such that for every interval I with |/| < 1, 

nii) < c\ir 

By an important theorem of Last |Las96| , UaH measures may be obtained 
by a process of closure. In our one-dimensional model the Hausdorff dimension 
varies locally in the s.c. spectrum, and the local Hausdorff dimension (suitably 
defined) may be determined explicitly. There exists a dense set in the s.c. 
subspace such that the spectral corresponding spectral measure is UaH. We 
shall for simplicity assume that the Hausdorff dimension has a constant value 
which will be denoted by a. 

The basic theorem, due to Strichartz |Str90] (with a very slick alternative 
proof by Last |Las96j . therefore we call it Strichartz-Last theorem), relates decay 
in the Cesaro sense to the Hausdorff dimension: 

Strichartz-Last theorem 

Let /i be a finite UaH measure and, for each / g £ 2 (R, dfi), denote 

fn(t) = J exp(-ixt)f(x)dfi(x) 

Then there exists C depending only on fi such that for all / G L 2 (R, cfyi) and 
T > 0, 

(|/ M | 2 ) T <C||/|| 2 T-" 

where ||/|| denotes the L 2 - norm of /, and (g)x = ^" 9 ^ dt . 
A measure is called a Rajchman measure iff 

lim fj,(t) = 

|t|->-oo 

It does not follow from the decay in the Strichartz-Last theorem that the corre- 
sponding measure is Rajchman . Let E denote the usual Cantor "middle-thirds" 
set in [— 7r,7r]; die F.S. transform of the corresponding measure T is: 

oo 

T(u) = IJcos[2/3u7t3^ +1 ] 

The corresponding Hausdorff dimension is well-known to be a = j }°| und 
the measure is UaH, but, from the above formula for T it follows that 

T(n) = r(3n) for all n e Z 

and hence the Cantor measure is not Rajchman. 
We are now ready to state our main result: 
Theorem 2 Let 

v 2 <o(v^-l) <v 2 c 
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with a < 4, where v c = 2-^/(/3 — 1). Then, for almost every (w 1 ,^ 2 ,^) with 
respect to v x v x I, [a.] there exist A with A + = — A~ and 

< A+ < A+ 

such that 

(a-(1 + 0),A+(1 + 0)) CfTac 

[b.] 

[-2(1 + 9), A-(l + 9)) U (A+(l + 0), 2(1 + 9)] C a pp 

[a] 

<Tsc(4 2) )n(\-(i + 9),\ + (i + 9)) 

may, or may not, be an empty set. 

The basic idea of the proof is that the Strichartz-Last theorem suggests a 
pointwise decay of the F.S. transform of the spectral measure of type t~ a ^ 2 for 
large t. For the Kronecker sum (with 9 = 1) the F.S. transform is the product of 
the corresponding transforms for the one-dimensional system, which we expect 
decays as t~ a . In order to use the proposition, we restrict the spectral measure to 
an interval in the s.c. spectrum such that 2a > 1 which is in principle non-empty 
(and may be proven so), and denote the F.S. transform of the thus restricted 
measure by the same symbol as before. It turns, however, out that this heuristics 
is not correct mathematically: as shown above, Cesaro decay does not in 
general imply pointwise decay or, in other words, that the spectral measure is 
Rajchman. This is a much harder problem - (see |MW12b] . chapter 4 for the 
full treatment of the pointwise decay of a model with superexponential sparsity 
and pure s.c. spectrum, and the recent very hard analysis to prove local decay 
in nonrelativistic QED TJFS12 ). Moreover, the result for general 6 may not be 
true due to resonance between Cantor sets, a subtle phenomenon which has been 
so far analysed only in the self-similar case by Hu and S.G. Taylor [HT94j . Self- 
similarity occurs, however, seldom: indeed, the s.c. spectrum of sparse Jacobi 
matrices is not self-similar by a theorem of Combes and Mantica CMOl]. We 
proved, however, that the main idea is indeed correct by generalizing a method 
due to Kahane and Salem ( jKS94j . |KS58j ) in specific (self-similar) cases to the 
problem at hand. 

Conceptual relevance of the present model: comparison with the 
Anderson model on a Cayley tree 

Sparse models in dimension d > 2 may be good models for the Anderson 
transition in lightly doped semiconductors ( which takes place already for d = 2 
): for the latter see [SE84]. It seems natural to expect that the present model 
might pave the way for a good qualitative description of the Anderson transition 
in lightly doped semiconductors. We say " pave the way" because the present 
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form of the model is not adequate for a physical description - but we argue that 
the main objection may be eliminated by considering a truly d -dimensional 
model. 

The main objection is, of course, that exponential sparsity is too severe, 
and not physically reasonable. It must be recalled, however, that the separa- 
ble model does not take account of dimensionality in a proper way, because 
"truly" three dimensional sparse models may drastically change, the cardinal- 
ity of {i : \a,i\ < R} in (4) from O (\ogR) to 0(R d ~ e ) in dimension d, for some 
e > 0, which is still compatible with the assumption of sparsity, changing, at 
the same time, the conditions on the sparsity for the existence of the transition. 

It should be useful to compare the present model with the Anderson model 
on a Cayley tree f |Kle98j . [JitOT] . [AWllj . |AW06Q . A Cayley tree (or Bethe lat- 
tice) is a graph with no closed loops: there is a central site and three generations 
in a regular tree of ramification r — 2. All sites are always connected to r + 1 
nearest neighbors. For a tree with N generations, there are N s = (r + l)r N ~ 1 
sites at the surface (generation N = 0) and a total Nt = (r N+1 + r N ~ 2 )/ (r — 1) 
sites. Since the ratio N$/Nt is nonzero for N — »• oo, the behaviour of statistical 
systems on a Cayley tree is usually (but see later) pathological and quite distinct 
from the physical features of a Bravais lattice. One notable example is the a.c. 
spectrum found in the Anderson model on a Bethe lattice in the seminal work 
of A. Klein [Kle98] (see also |Jit07j for a review and further references): the 
extended states decay exponentially but are not square-integrable due to the 
exp(n) many points within n links of a given vertex alluded to before, and is 
thus of entirely different nature from the a.c. spectrum found on a Bravais 
lattice - and in particular in theorem 2. 

It is often mentioned in the literature that the Bethe lattice describes the 
infinite-dimensional limit of a hypercubic lattice. Indeed, the Bethe or Bethe- 
Peierls approximation (BPA) ([Hua87], pg. 357), which improves on the mean 
field approximation (which describes the infinite-dimensional limit) by taking 
into account specific short-range order, agrees with the Bethe lattice in the fol- 
lowing sense |Tho82j : expectation values of spins "far removed from the surface" 
on a Cayley tree agree with those obtained by the BPA. That this is truly only a 
local property is made clear by the fact that, by first considering a finite Bethe 
lattice and then performing the thermodynamic limit one does not obtain the 
BPA (see again the remarks in )Tho82j ). 

Klein's proof of the Anderson transition |Kle98| uses the loopless character of 
the graph in an essential way and thus, as remarked by Jitomirskaya [Jit07] , the 
Bethe lattice, while infinite-dimensional (but in the above-mentioned sense), is, 
in a sense, one-dimensional. Similarly, our model has no loops, because theorem 
2 "inherits" the one-dimensional structure of the model in theorem 1. The 
geometrical structure underlying (6), i.e., the three coordinate half-axes (there 
would be no problem in extending it to cover the three full axes) is, however, 
a (small!) part of the full (sparse) model on a Bravais lattice. Thus, in our 
opinion, in spite of the considerable independent mathematical interest of the 
recent work on the Anderson model on a Bethe lattice ( |AW11| . jAW06j ) . the 
present model is physically more reasonable (for the description of the Anderson 



10 



transition in lightly doped semiconductors). Of course, the flaws discussed at 
the beginning of this subsection will only be disposed of by "filling in" the 
"remaining" points -an immensely challenging, basic open problem!. 

As in the Bethe lattice [Klc98], the sharpness of the transition, i.e., the 
existence of a mobility edge, was not proved for the present model. Recent 
work [AWllj proves that no mobility edge occurs in the Bethe lattice at weak 
disorder. The general character of the arguments used in theorem 2 to establish 
the existence of a.c. spectrum (see the sketch of the ideas of the proof there and 
|MW12aj ) suggests that the intermediate region might be more accessible to 
analysis than the Bethe lattice, but this remains a (challenging) open problem. 

It is, however, rewarding that already the separable model displays a dra- 
matic "kinematic" effect of the dimensionality: for d > 2 the transition becomes 
truly Anderson-like, i.e., from a.c. to p.p. spectrum, we say dramatic because, 
in the case of heavy doping, the d > 2 version, built as in (6) from the one- 
dimensional version of (1) - (3), continues to have purely p.p. spectrum, by 
[GMP77] and the same proof of theorem 2 (more precisely, the easy part of the- 
orem 2), in complete disagreement with the expected transition! Thus, theorem 
2 is a definite indication that the present approach via "light doping" is more 
likely to produce a transition in the full (sparse) version. 

Frustration and short-range spin glasses: the Edwards- Anderson 
model 

Dilute solutions of atoms of large magnetic moments (such as the transition 
metals Fe, Co, Mn) in a paramagnetic substrate (Cu, Au) present a number 
of peculiar physical properties. For small, but sufficiently high concentrations 
of magnetic impurities, the susceptibility in low fields displays a characteristic 
peak, with discontinuous derivatives, at a temperature T g . The specific heat is 
always smooth, with a linear dependence on the temperature as T — > 0. The be- 
haviour of these spin glasses has been explained in terms of an indirect RKKY 
interaction between the spins of these magnetic impurities mediated by the elec- 
trons of the paramegnetic matrix [BY861 . which is of the form — J(\i — j\)SiSj 
with J(r) = (kpr)~ 3 cos(2kpr), where r is the distance between magnetic atoms 
and kp the fermi momentum; Si = ±1 are (e.g.) Ising spins. The rapid oscil- 
lations and the weak decay of J(r), as well as the random distribution of the 
impurity magnetic atoms, are the basic ingredients of spin glasses. At sufficiently 
low temperatures there is a "freezing" of the magnetic moments in random di- 
rections (which leads to an increase of the susceptibility). The spin glass phase 
may be regarded as a conglomerate of blocks of spins, each block with its own 
characteristic orientation, in such a way that there is no macroscopic magnetic 
moment [BY86]. 

Edwards and Anderson |EA75j proposed a spin Hamiltonian, which we write 
in a generalized version as follows. For each finite set of points A c Z d , where 
the dimension d will be restricted to the values d = 2 and d = 3, consider the 
Hamiltonian 




(7a) 



ijeA 
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where <&ij with i, j £ A are self-adjoint elements of the algebra generated by 
the set of spin operators, the Pauli matrices erf, a?, erf, i € A, on the Hilbert 
space "Ha = ®igACf, given by 

$*J = + a a cr,f crj + <\:<t',(t] (7 - 6) 

for \i — j| = 1, and zero otherwise. The random couplings Jij, with \i — j\ = 1, 
are random variables (r.v.) on a probability space (Q, J 7 , P), where J 7 is a sigma 
algebra on f2 and P is a probability measure on (fi, J 7 ). We may take without 
loss of generality 

n=x Bd S (8) 

where S 1 is a Borel subset of R, B d is the set of bonds in d dimensions, and 
assume that the Jjj are independent, identically distributed r.v.. In this case, 
P is the product measure 

dP=x Bd dP (9) 

of the common distribution Pq of the random variables, which will be denoted 
collectively by J. The corresponding expectation (integral with respect to P) 
will be denoted by the symbol Av. We have to assume that 

Av(J id ) = (10) 

for all i,j £ Z d , i.e., that the couplings are centered. This assumption mimicks 
the rapid oscillations of the RKKY interaction. Let Ea denote the GS energy 
of H\, i.e., E\ = inf spec(HA). The following result was proved, among several 
others, in |CL10j : 

Theorem 3 ( [CLIO] ) For P- a.e. {J}, the limit below exists and is inde- 
pendent of the b.c: 

e <«0= ii m (ii _ a ) 

A/Z" |A| ' 

and 

e ( rf )=mf^ (11-6) 
a |A| 

where |A| denotes the number of sites in A. Finally 

e W > e ( d+1 ) (11 - c) 

(11-a) is the far-reaching property of self-averaging (see .An d78] for a dis- 
cussion): it expresses that measurable - e.g. thermodynamic- quantities are 
the same for any typical configuration of the sample, i.e., are experimentally 
reproducible. It follows from (11a) that, P- a.e., 

eW = lim Av &*) (11 _ d) 

A/*Z* |A| ' 

Let Aat denote a square with N sites if d — 2 or a cube with sites if d = 3, 
and write H^\j) = H\ N (J). We now adopt periodic b.c. for simplicity, but 
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the final result is independent of the b.c. due to theorem 1. We may write 



where is given by 



H^{J)=c d V JtjQij (12 -b) 



Above, Cd are factors which eliminate the multiple counting of bonds, i.e., 

c 2 = 1/2 (12 

and 



c 3 = 1/4 (12 - d) 

and A„ is a square labelled by a site n, for which we adopt the convention, 
using a right-handed (x, y) coordinate system, that n = (n x ,n y ) is the vertex 

in the square with the smallest values of n x and n y . Similarly, A„ is a cube 
labelled by a site n = (n x ,n y ,n z ) with, by the same convention, the smallest 
values of n Xl n y and n z . Due to the periodic b.c, the sum in (12-b) contains 
precisely N lattice sites. The notation is short-hand for its tensor product 
with the identity at the complementary lattice sites in An\A^\ Let EQ N ' d \j) 
denote the GS energy of H$ (J) , E [ ™ 4) {J) the GS energy of i#"(J) and 

E^ d) = Av{E { *> d \j)) (12 -e) 

By the condition of identical distribution of the r.v. Jy, E^ does not depend 
on n, which is implicit in the notation used. We have 

Theorem 4 (Theorem 2 and proposition 1 of |Wrel2j ) The following lower 
bound holds: 

MEf«M iS p (13 _ 0) 

Further, let dPo in (9) be the Bernoulli distribution dPo = 1/2(8 j + 5-j), set 
for simplicity J = 1 and let a x — a y — in (7b). Then, 

> -3/2 (13-6) 

r ,i , ,36096 

e (3) > _x/4 (13 -c 

- 1 4096 v ' 

The special case a x — a y = in (7-b) is the classical EA spin-glass; we also 
set a z = 1. In this case e^ d \ given by the r.h.s. of (11-d), is invariant under the 
"local gauge transformation" erf — > —erf together with J,j — > — Jij for all j\\j — 
i\ = 1, whatever the lattice site i. 



and 
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According to the above, |Tou77j . an elementary square (" plaquette" ) P is 
said to be frustrated (resp. non-frustrated) if Gp = j)ep^ij = — 1 resp. 
+1. Note that Gp is gauge-invariant and that, for the quantum XY (or XZ) 
model defined by setting a y — in (7-b), is also locally gauge-invariant if 
we add to the above definition the transformation erf — » —erf, i.e., the trans- 
formation on the spins is defined to be a rotation of tt around the y - axis in 
spin space. The property of local gauge-invariance guarantees the absence of a 
macroscopic magnetic moment (spontaneous magnetization) mentioned before 
as a basic property of spin glasses [J ASM] ■ 

Since the Pauli z-matrices commute, finding the minimum eigenvalue of (7- 
a) in the classical EA case is equivalent to find the configuration of Ising spins 
<r, = ±1, denoted collectively by a, which minimizes the functional 

F(a, J) = Jijcno-j 

a 

The minimal energy of a frustrated plaquette equals Epj — —2 and of a non- 
frustrated plaquette Ep. u = —4. 

(13b, c) of theorem 4 provide the first (nontrivial) rigorous lower bound both 
for e^ 2 ) and e^ 3 ^. Using the natural misfit parameter 

\m - \ En \ 



\E*\ 

as a measure of plaquettes frustrated or bonds unsatisfied (see (4) of [KK95J, 
where Eo denotes the ground state energy of the frustrated system and E ld is 
the ground state energy of a relevant unfrustrated reference system, we find from 
(13b) in the d = 2 case the lower bound m > 0.25 and for d = 3, from (13c), 
the lower bound m > 0.26 • • • : thus, in both cases, the measure of frustrated 
plaquettes or unsatisfied bonds as defined above is at least of the order of 1/4. 

The method of proof of theorem 4 is a rigorous version of a finite-size cluster 
method, originally due to |BS79j . together with the variational principle. If we 
take in (7) a z = 1, a y = 0, and a x = a, and consider a as a small parameter, 
we have the anisotropic XZ (or XY) model. By the norm-equicontinuity (in the 
volume) of H\({J}) (given by (7)) as a function of a, which is preserved upon 
taking averages over the probability distribution of the J, it follows that (13b, c) 
hold with the right hand sides varying by small ammounts if \a\ is sufficiently 
small: this is conceptually important for reasons mentioned in the introduction. 

The mean field theory, recently rigorously solved in a spectacular way (see 
|Gue06] for a review) does not exhibit frustration - indeed, this concept is not 
even generally defined for their model, since the theory does not require an 
underlying lattice. Whether frustration is an important issue in the description 
of realistic spin glasses is an important open problem. 

We refer to the article by Bovier and Frohlich [BF86 (see also Bovier's book 
|Bov06| ) for an illuminating discussion of complementary, mostly global (i.e., 
involving the lattice as a whole) aspects of frustration. Our bound (13c), which 
relies only on the local structure, is, however, slightly better than Kirkpatrick's 
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|Kir77j . which is based on very reasonable, but unproved conjectures of a global 
nature. A most relevant open problem would be, of course, to extend the finite- 
size cluster method to obtain bounds for the free energy of the EA model. 

Return to equilibrium and (thermodynamic) instability of the first 
and second kinds 

The dynamics of spin glasses is a topic of great relevance, both conceptual 
and experimental. It happens, however, that the standard approach to the 
kinetic theory (e.g., for the mean-field spin glass model) relies on Glauber dy- 
namics (see, e.g., |Sza97] ).a well-known dynamics imposed on the Ising model, 
which has been used to study metastability in early days [DC 074] and more 
recently [SS98J. In spite of the considerable independent mathematical interest 
of the latter works, it should be remarked that only for a quantum system does 
a physically satisfactory definition of the dynamics of the states and observ- 
ables exist which is relevant to the microscopic domain, including, of course, 
condensed matter physics. The scarcity of examples of approach to equilibrium 
in quantum mechanics is due, of course, to the extreme difficulty of estimating 
specific properties of the quantum evolution. 

This fact was the basic motivation which prompted Emch [Emc66] and Radin 
[RadTQ] to propose a quantum dynamical model (which we dubbed the Emch- 
Radin model in |Wrel2j and MW12b ), which is relevant to systems with high 
anisotropy, and displays a remarkable property of non-Markovian approach 
to equilibrium, or return to equilibrium. It turns out that such a model is 
useful not only for a description of nuclear spin-resonance experiments - such 
as the one |LN57| which motivated Emch and is described below - but also to 
describe dynamical effects associated to quantum crossover phenomena in spin 
glasses (see |Sac94| for a review and references). For this purpose, a material 
was chosen with a strong spin-orbit coupling between the spins (of the magnetic 
ion) and the underlying crystal: this coupling essentially restricts the spins to 
orient either parallel or antiparallel to a specific crystalline axis, which we shall 
label as the z-axis. Such spins are usually referred to as Ising spins, and will be 
described by the Ising part of the forthcoming Hamiltonian. In the experiments 
(see Sac94j and references), a transverse magnetic field is then applied, oriented 
perpendicular to the z axis. A large enough transverse field will eventually 
destroy the spin glass order even at T = 0, leding to the existence of a crossover 
region. 

We shall be interested in a situation in which a large transverse field has 
been applied; the initial state of the system is, then, approximately, a product 
state of the forthcoming form (16). The same experimental setup should there- 
fore allow a measurement of the rate of return to equilibrium of the mean 
transverse magnetization, i.e., of how f(t), defined by (20), tends to zero, or 
of how (21) is approached. As we shall see, this rate depends strongly on the 
probability distribution of the J. 

The XY model has also been extensively studied from the point of view of 
return to equilibrium, since early days as an example of general theory |Rob73j , 
more recently in AB83] (see the latter for reference to related important work 
of Araki) , but unfortunately only in one dimension. 
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The nonrandom model 

Following Radin's review [Rad70| of the work of Emch |Emc 66j . consider 
the experiment of [LN57 : a CaF 2 crystal is placed in a magnetic field thus 
determinig the z-direction, and allowed to reach thermal equilibrium. A rf pulse 
is then applied which turns the net nuclear magnetization to the x direction. 
The magnetization in the x direction is then measured as a function of time. 

As in f[Emc66 . Rad70]h we assume an interaction of the form 

H V = 1/2^^3 -k\)a*ai (15) 

where j € V and k € V in (15), with V is a finite region, V E Z . Hy is defined 
on the Hilbert space "H — C 2 . The state representing the system after 

the application of the rf pulse will be assumed to be the product state 

P = P(0) = (g) fa (16 - a) 

where 



^■(•) = trj(- exp(- 7 aj ; ))/tr J (exp(-7^) (16 - b) 

and tr is the trace on C 2 . Other choices of the state are possible |Rad70j , but 
we shall adopt (16) for definiteness. Let 

S* = 1/N(V) J2 (17) 
jev 

be the mean transverse magnetization, with N(V) denoting the number of sites 
in V. The real number 7 in may be chosen as in |Emc66] such as to maximize 
the microscopic entropy subject to the constraint trS x p(0) equal to a constant, 
i.e., a given value of the mean transverse magnetization. Since the state (16) is 
a product state, 7 is independent of V , and has the same value if S x is replaced 
by af , for any i € V. Let 



and 



where Uv(t) = exp(iiiTy). We may write, by (18) 



p t =WO)¥ (18- 



<>y(t) = pr(<) (18-6) 



)v(t) = p(0)(U v (-t)a? o U v (t)) (18- 



It is natural to define 



«>(*) = lim (af B ) v (t) (19) 

provided the limit on the r.h.s. of (19) exists, as the expectation value of the local 
transverse spin in the (time-dependent) nonequilibrium state of the infinite 
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system. A weak-star limit of the sequence of states p v (-) = tr(-U v p^Uy 1 ) 
exists, by compactness, on the usual quasilocal algebra A of observables (see, e.g. 
|Sew86j . for a comprehensive introduction to the algebraic concepts employed 
here). The expectation value of af o in the equilibrium state associated to (15) 
is zero by the symmetry of rotation by it around the z axis. We may now pose 
the question whether the limit 

/(t)= lim pL(lMA)$>f o ) (20) 

where A denotes a finite subset of Z", which is interpreted as the mean transverse 
magnetization, exists. The property of return to equilibrium is expressed by 

lim f(t) = (21) 

t— >oo 

Of particular interest is the rate of return to equilibrium. In the present non- 
random case, the limit at the r.h.s. of (21) equals /(<) = /^(cf ), for any io, 
by translation invariance of p^. This is not so for random systems, in which 
case additional arguments are necessary to show the convergence of the r.h.s. 
of (21) for almost all configurations of couplings. 
We assume that the function e(-) satisfies: 

e(0) = (22) 

and either: 



^ e(j) < oo i.e. e e Zi(Z d ) (23 - a) 

]£Z d 



or 



£ (e(j)) 2 < oo i.e. 6 G h(Z d ) (23-6) 
The random model 

We now introduce the Hamiltonian of a disordered system corresponding to 
(15), where we take, for simplicity, 

V = V n = [-n,n} d (23 -d) 



fl» = l/2 J{iMl- k K°l ( 24 ) 

j,kev n 

where J{j,k) arc independent, identically distributed random variables (i.i.d. 
r.v.) on a probability space which we denote by P). We shall use Av(-) 

to denote averaging with respect to the random configuration { J^fe}, which we 
denote collectively by the symbol J. The J are assumed to satisfy |vEv H83j : 

Av(J m ) = Q (25 -a) 

\Av(J^ k) )\<n\c n Vn = 2,3A,--- (25-6) 
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We shall take, without loss, e(j) > 0. If 
<J) - 



ifjeifc, 

otherwise , 



(25 



where ±8i for i = !,■■■ ,d denote the set of z = 2d bonds connecting the origin 
to a point of Z d , we have the EA spin glass model of section 3. 

Let \V n \ denote the number of sites in V n and the free energy per site f n be 
defined by 

~kT\ogZ n (J) 
fn(J) = |^rj (26 - a) 

where 

Z n (J)=tr(exp(-PR n )) (26-6) 

is the partition function and the trace is over the Hilbert space T~L. Then 

Theorem 5[vEvH83 (see also [KS79 for the first result in this direction): 
Under assumptions (12-b) and (25), the thermodynamic limit of the free 

energy per site 

f(J) = lim f n (J) (27 -a) 

n— too 

exists and equals its average: 

f(J) = Av(f(J)) = lim Av(f n (J)) (27 - b) 

n— > oo 

for almost all configurations J ( a.e. J). 

The reason why (23-b) suffices for the existence of the thermodynamic limit 
is that, in order to obtain a uniform lower bound for the average free energy 
per site, the cumulant expansion (see, e.g., ( Si m79| . (12.14), pg 129, for the 
definition of Av c , there called Ursell functions): 



n=2 

was used |vEvH83] . which, by (25-a), starts with the second cumulant 
A« c (JEj-0 = Av{ j£ jj), which is the variance of J(ij). Condition (25-b) was 
used to control the sum in the exponent above. 
Thermodynamic stability 

Stability considerations play a key role in quantum mechanics [LSlOj . Let 
a statistical mechanical system be described by a collection of amiltonians 
associated to finite regions A C R d , for particle systems, or A C Z d for spin 
systems, with volume V(A): examples are (15) and, for random systems, (24): 
in the latter case there is implicit in the Hamiltonians a dependence on the 
random variables J, and the constant c in definition 2 below is assumed to be 
a.e. independent of J. The system's free energy is /a = ~ fcT ^ s A ^ A< -' 7 - ) with 
Z\ = tr (exp(— /? Ha)), and the thermodynamic limit means A /*• R d (or Z d ) 
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with the proviso of fixed density for particle systems, where A /■ Z d denotes a 
limit in the sense of van Hove or Fisher (the latter being required for random 
systems, see the discussion in |vEvH83] : (23-d) satisfies this assumption). 

Definition 2 The system is said to be thermodynamically stable if there 
exists a constant < c < oo such that 

/a > -c (28 - a) 

It is said to satisfy thermodynamic stability of the first kind if 

H A > -oo (28 - b) 

and to satisfy thermodynamic stability of the second kind if 

H A > -cV(A) (28 - c) 

The above definitions (28-b,c) are patterned after the corresponding ones for 
N- body systems in |LS10| : they also apply to relativistic quantum field theory, 
where the particle number N is not conserved. 

Theorem 5 illustrates the interesting fact that, for disordered systems, sta- 
bility of the second kind (28-c) may fail even for a thermodynamically stable 
system: this happens when only (23-b) (but not (23-a)) holds. In the following, 
we shall demonstrate that this phenomenon has important dynamical impli- 
cations for a relevant class of disordered systems. 

We now return to our nonrandom model (15). When (23-a) holds, it follows 
from a folklore theorem (see, e.g., Theorem 6-1 of [MW12b] or Theorem 3.3, pg. 
Ill, of |Dav76| ) and a representation of / ( [Rad70j . pg.295, and proposition 1) 
that exponential decay in the sense that 

\f(t)\<Cexp(-d\t\) (29) 

with C and d positive constants, cannot hold. This is essentially the condition 
that the physical (GNS) Hamiltonian H is bounded from below (semibounded) 
(see again |Sew 86 ). where 

Ho = <r z o < k ) a k ( 30 - °) 

kez d 

The infinite sum (30) stands for a limit in norm in the quasi-local algebra A 
associated to the spin system, see Rad70 , loc. cit. It is the thermodynamic 
stability of the second kind (28-c) which guarantees the existence of Hq, and 
(29) as a consequence. In the random case (24), 

Ho(J) = <?o J o,fc £ ( fc K (3° - b ) 

kez d 

under assumption (23-a): if only (23-b) holds, the representation (30-b) is, of 
course, not defined. We refer to |Wrel2] . proposition 4.1 for the dynamical 



19 



consequences of this fact: in particular, unlike the stable case (23-a), exponen- 
tial decay (29) does indeed hold for a class of one-dimensional potentials with 
algebraic decay, as shown there. 

We now consider the special nearest-neighbor interaction (25-c) - i.e., the 
EA model of section 3. The following theorem may be proved (along the same 
lines, but simpler, than theorem 3.2 of |Wrel 2 )[^| 

Theorem 6 Let (25-c) hold, i.e., (24) describe the EA spin glass model 
of section 3. Then there exists a subsequence {m n }nez such that, a.e. with 
respect to J, 



f(t)= lim pU 



where 
and 



= SAv(p (t)) 

(31) 

S = 0oK) (32) 
p (t) = []cos(2/3 J ,± Si t) (33) 

±5, 

We have the 

Corollary 6 Consider the distribution functions 

dP^x) = l/2[<5i(a;) + 8-i{x)] Bernoulli distribution (34 - a) 

dPijdx = l/2x[_ 11 ](a;) uniform distribution (34 — b) 

dP^/dx = — -j= — cxp(— x 2 ) Gaussian of unit variance (34 — c) 

The corresponding values of /, defined by the r.h.s. of (31), are: 

fx{t) = (cos(2f3t) z (35 -a) 

Ht) = {—^- L r (35-6) 
/ 3 (t) =cxp(-2z< 2 ) (35 -c) 

Remark 2 

a. ) In the Bernoulli case, (35-a) implies no decay, as in the non-random 
model: this is an instance of a result of Radin for general interactions satisfying 
(23-a), according to which / is almost-periodic ([Rad7(J, Lemma 2, pg.1951). 

b. ) For the uniform distribution, (35-b) describes a non-Markovian re- 
turn to equilibrium. Clearly, in this case, under assumption (23-a), (30b) is 



1 There are two misprints in the latter reference: in the r.h.s. of (40a) the number four 
should be inside the exp function, while in (40b) the first parenthesis should be moved to the 
left of the symbol Av, and the last parenthesis omitted. 
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a physical Hamiltonian which is semiboundcd for each realization of J, and thus 
exponential decay (29) cannot hold by the aforementioned folklore theorem. 

c. ) In the Gaussian case, not even thermodynamic stability of the first kind 
(28-b) holds, because the Gaussian r.v. range over R. (30b) shows that also the 
physical Hamiltonian H is not bounded below, and thus (29) may, in principle, 
hold. That this is actually the case is, of course, made clear by (35-c). 

d. ) Since the transverse magnetization in (31) is an observable (measurable 
quantity), self-averaging is an important requirement, as discussed before. We 
conjecture that all subsequential limits are equal, a.e. J, to the r.h.s. of (31), 
i.e., that the limit n — > oo in (31) exists. 

The issue of probability distributions 

As is well-known, the use of Gaussian probability distributions (p.d.) in dis- 
ordered systems is standard: it simplifies several passages considerably, e.g., the 
integration by parts formula in the first proof of the existence of the thermody- 
namic limit for the mean field theory in |GT02| , or in the first proof of diffusion 
in the full Anderson model |FS83j (the proof here allowed also the uniform dis- 
tribution (34-b), considered by Anderson |An d58 ). One might think that the 
results are expected to depend qualitatively on the probability distribution: 
this is, however, not so. 

In |S W85j . where a mean-field Ising model in a random external field was 
studied, it was proved that the existence or not of a tricritical point in the 
model's phase diagram is tied to the probability distribution: it is absent for a 
Gaussian p.d., and present for a Bernoulli distribution. 

As an attempt to explain this result physically, we conjectured that it was 
due to the fact that a discrete distribution of probabilities such as the Bernoulli 
distribution samples just a few values of the couplings and thus introduces some 
short-ranged elements into the problem: from this point of view, discrete dis- 
tributions may have a closer connection with real materials. This feature is 
shared by the uniform distribution, see remark 2 b.), but even here there is a 
remarkable dynamical difference between Bernoulli and uniform distributions 
with regard to decay, exhibited by (35-a) and (35-b). The latter yields, some- 
what surprisingly, the same result in one dimension (z = 2) as obtained by 
Emch |Emc66| with an interaction of infinite range (e(|n|) = 2~l"l~ 1 ) in the 
nonrandom model (or, alternatively, the random model with Bernoulli distri- 
bution and the same potential). Since the latter (algebraic slow decay with 
oscillations) seemed to lead to a good qualitative description of the nuclear spin 
resonance experiment ILN57| . the uniform distribution decay (35-b) might well 
be in qualitative agreement with the previously proposed spin glass experiment! 

In contrast to a discrete or uniform distribution, a Gaussian distribution 
samples many values of the couplings and works to reinforce the long-range 
nature of the interactions. For effectively short-range inteactions, this might 
not be suitable. In any case, the Gaussian probability distribution leads to a 
spectacularly fast rate of return to equilibrium (35-c), and such sharp differences, 
such as between (35-c) and (35-b), should be amenable to experiment. 
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